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Abstract 


There are a number of existing methods to transfer information among various disciplines. For a 
multidisciplinary application with n disciplines, the traditional methods may be required to model (n 2 — n) 
interactions. This paper presents a unified three-dimensional approach that reduces the number of 
interactions from ( n 2 -n ) to 2 n by using a computer-aided design model. The proposed modeling 

approach unifies the interactions among various disciplines. The approach is independent of specific 
discipline implementation, and a number of existing methods can be reformulated in the context of the 
proposed unified approach. This paper provides an overview of the proposed unified approach and 
reformulations for two existing methods. The unified approach is specially tailored for application 
environments where the geometry is created and managed through a computer-aided design system. 


Results are presented for a blended-wing body and a 

Introduction 

A key element in the application of 
multidisciplinary design optimization (MDO) to 
an engineering system is the introduction of a 
consistent geometric representation. Such a 
representation guarantees that the same 
geometry model is used to derive the 
computational models required for various 
disciplinary analyses. By utilizing computer- 
aided design (CAD) for consistent geometry 
representation, it is easier to analyze complex 
configurations with higher-fidelity tools such as 
computational fluid dynamics (CFD), 
computational structural mechanics (CSM), or 
detailed finite-element analysis. 

A feature that characterizes multidisciplinary 
analysis and optimization is the modeling of 
interactions among various disciplines. For 


ugh-speed civil transport. 

example, the strong interactions between CSM 
and CFD can prompt physically important 
phenomena such as those occurring in aircraft 
due to aeroelasticity. Correct modeling of these 
complex aeroelastic phenomena requires a 
coupling of CSM and CFD for a flexible structure 
(e.g., airplane). In a multidisciplinary 

environment, various disciplines must represent 
the same configuration geometry, and data from 
each discipline must be available consistently to 
all the disciplines. The data may be scalar (e.g., 
pressure and temperature), vector (e.g., 
deflection and heat transfer), or integrated 
quantities (e.g., aerodynamic and thermal 
loads). The data transfer process may be 
subjected to additional constraints, such as 
conservation of forces, moments, and energy. 
The focus of this paper is the transfer of data 
between dissimilar grids (most models do not 
share the same nodal locations at the interface). 
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Review of Existing Methods 

The issue of exchange of pressure distribution 
and displacements between an aerodynamic 
model and a structural model gained urgency 
with the first widespread use of panel 
aerodynamic methods and finite-element 
structural models during the late sixties. The 
traditional beam and strip theory methods were 
not applicable to the large, flexible low-aspect 
ratio configuration of a supersonic transport 
(SST). The SST experience led to development 
of a comprehensive aeroelastic computer 
program called FLEXSTAB. 1 Most of the work 
was done under NASA funding before 1974. The 
basic theoretical work was completed between 
1968-71, and the goal of FLEXSTAB was to 
provide an aeroelastic program for production 
use. FLEXSTAB probably was the first 
computing system to address systematically all 
the aeroelastic transfer issues for a complete 
airplane configuration. FLEXSTAB was 
successful in providing an acceptable method 
for the pressure and displacement transfers. 
However, the method was time consuming, and 
required significant manual effort for the user to 
establish direct correspondence between the 
aerodynamic panels and structural elements 

In the eighties Dassault Aviation developed a 
more simplified method 2,3 including many 
practical considerations. The Dassault method is 
embodied in Elfini and includes several 
innovative ideas in aeroelastic formulation. First, 
Dassault used an intermediate computational 
grid between an aerodynamic mesh and a 
structural finite-element model for a given 
configuration. Thus different aerodynamic 
models could be used with the same structural 
model, and vice versa. Second, they used shape 
functions defined on the computational grid, and 
smoothing operators to relate structural 
displacements to the aerodynamic mesh. Third, 
they introduced the idea of loads basis on the 
computational grid (unit loads at the nodes of 
the computational grids) to transfer pressure 
from the aerodynamic mesh to the 
computational grid, and from the computational 
grid to the finite-element nodes. The overall 
approach is elegant and practical. But it still 
requires significant user input (although much 
less than FLEXSTAB) and checking to ensure 
the accuracy of the transfers. 


In the recent years, various researchers have 
examined the issue of aeroelastic transfer. 
Discrepancies and dissimilarities in geometry 
and grid models are two potential sources of 
error. The accuracy of the data transfer depends 
on the relative resolutions of disciplinary grids. 
Data could be lost in transfer from a coarse grid 
to a fine grid. Another source of error occurs if 
the models have dissimilar levels of geometry 
detail. As noted by Tzong et al., 4 a CFD grid 
generally resembles the true geometry of the 
aircraft; the grid includes details such as pylons, 
nacelles, flaps, and slats. However, a CSM grid 
generally represents only major structural 
components, such as the wing box. Flaps and 
slats are represented either by a few simple 
beam elements or are completely excluded. 
Tzong et al. 4 and Kapania and Bhardwaj 5 have 
developed methods, based on finite element 
(FE) technology, in which virtual work is 
employed to transfer the aerodynamic pressures 
onto a CSM grid. The displacements are then 
converted back to a CFD grid through the 
reciprocal theorem. Kapania and Bhardwaj 5 
were successful in using a simplified version of 
this method for several wings. 

Brown 6 added virtual elements in the CSM 
model to cover the discrepancies in the 
geometry definition between CSM and CFD 
models. These virtual elements add neither 
stiffness nor mass. As pointed out by Cebral and 
Lohner, 7 the generation of the virtual elements 
is an unnecessary complication, particularly for 
complex geometries. 

Hounjet and Meijer 8 and Smith, Hodges, and 
Cesnik 9 provided overviews of the data transfer 
methods. Smith, Hodges, and Cesnik 9 evaluated 
six methods for transferring information between 
CFD and CSM disciplines. These methods were 
infinite-plate spline (IPS), multiquadric 
biharmonic (MQ), nonuniform B-spline (NUBS), 
thin-plate spline (TPS), finite-plate spline (FPS), 
and inverse isoparametric mapping (IIM). These 
methods have been implemented in a single 
code, FASIT 10 . The IPS method is based on the 
popular surface splines 11 and is implemented in 
some commercial aeroelastic analysis tools. 
This method is designed for interpolating a 
function of two variables. 

Out of these six methods, Smith, Hodges, and 
Cesnik 9 recommended further study of IIM and 
NUBS. They indicated that IIM shows great 
promise for two-dimensional applications and 
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needed to be extended to three dimensions. 
Clutter' 2 and Send 13 extended NUBS to three 
dimensions. 

One problem found with the NUBS 
implementation is that the data must be input as 
a structured (regular) grid. This requirement 
forces the data, at best, to be approximated, and 
in most realistic cases, this step is either time 
consuming or impossible. Samareh 14 proposed 
a method to use non-uniform rational B-spline 
(NURBS) representation for data transfer among 
various disciplines. Because this method is 
based on a general three-dimensional, least- 
squares representation, 15-16 it does not require 
the input to be a structured grid. Another 
advantage of this approach is control over the 
tradeoff between smoothness and accuracy. 

The accuracy of the data transfer process for 
integrated quantities (e.g., forces, moments, and 
energy) depends on the consistency of data 
transfer as well as other constraints, such as 
conservation. For example, a consistent load 
vector for CSM is defined as 


proposed by Murti 19 for a two-dimensional 
model; it uses the FE shape functions to 
interpolate the coordinates, pressure, and 
displacement vectors. The method does not 
require any matrix inversion. Because the FE 
shape functions satisfy a positivity constraint, 
the process will not create nonexisting local 
extrema. The local interpolation is computed by 
projecting one grid onto another. Cebral and 
Lohner 7 ' 26 presented a variation of the I BA that 
could guarantee the conservation of forces. 
They used a Galerkin method to solve for 
pressure from the CSM grid. The Cebral-Lohner 
method requires a matrix inversion. They also 
used an adaptive Gaussian integration 
technique to improve the accuracy. Farhat, 
Lesoinne, and LeTallec 21 also presented a 
variation of the original IBA 18 that can guarantee 
conservation of forces, but requires no matrix 
inversion. We used this algorithm for the current 
study. 


Overview 


{/ e } = J{W} r n{FWS 

s 

where {f e } is the element load vector, {A} is 
the FE shape function, n is the unit surface 
normal, P is the pressure, and dS is the 
infinitesimal surface element. 


The unified approach has two essential 
ingredients. First, the data transfer process 
between two disciplines was modeled by a 
transformation matrix. Second, the CAD model 
was used to reduce the number of interactions. 

The interaction between two disciplines is 
modeled mathematically as 


Because the above equation uses the same 
shape functions as are used to calculate the 
element stiffness matrix, the equation 
guarantees a consistent loading or lumping. 
Cook, Malkus, and Plesha 17 provided a detailed 
description of this equation. There are several 
possible problems with using the above 
equation. First, the aerodynamic load may have 
a large variation within a single FE, such that the 
shape function is not adequate to capture the 
variation. Second, the FE shape function may 
not be available for some commercial CSM 
codes. Third, in its present form the above 
equation does not guarantee conservation of 
forces and moments. 

Maman and Farhat 18 outlined a consistent 
interpolation-based algorithm (IBA) for 
transferring information between two dissimilar 
grids. The algorithm is similar to the IIM 


{F 2 } = [T 2 l ]{F ] } 

where matrices {F,} and {F 2 } contain the 
information on discipline grids 1 and 2, 
respectively, and matrix IT 2I ] is a 
transformation matrix. For example, {F,} could 
be an aerodynamic loads vector defined on a 
CFD grid and transferred to the CSM grid as 
{ F 2 } . Generally the transformation matrices 
are sparse and large. Only the nonzero 
elements need to be stored. If the transformation 
matrix [F 21 ] is independent of the shape 

changes, then [F 21 ] can be calculated once 

and used as long as there is no change in the 
grid connectivity. 

The concept of a transformation matrix simplifies 
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integrated analyses such as aeroelastic 
calculation. The aeroelastic calculation has four 
distinct steps. First the aerodynamic loads are 
calculated on the CFD grid. Second the loads 
are transferred to the CSM grid. Third the 
aeroelastic deflections are calculated on the 
CSM grid. Fourth the deflections are transferred 
to the aerodynamic grid to recalculate the 
aerodynamic loads. This iterative process can 
be expressed as 

{f f }= {Flow Solution ( G F + <5 F )} 

Mfc}= ft) 

The first equation represents the aerodynamic 
load calculation. The term G F represents the 

CFD grid, and S F is the aeroelastic deflection 
on the CFD grid. The aerodynamic forces, F f , 
are transferred to the CSM grid as F s . The 
matrix [F] is the CSM stiffness matrix, and 
S s is the aeroelastic deflection on the CSM grid. 

The use of the transformation matrix simplifies 
the above set of equations to 

fc}=fc(c F + [r ra tf,5 s })} 

For a linear structure without rigid body degrees- 
of-freedom, the above equation set can be 
simplified to 

where 

[r FF Hr ra M'[rJ 

The formulation can be extended to the case 
with rigid body degrees-of-freedom, and it allows 
decoupling of the CFD code from a linear CSM 
code. 

If the transformation matrix is independent of the 
shape changes, then the formulation is 
especially beneficial. The benefit results from 


the fact that the sensitivity derivatives used in 
gradient-based optimization with respect to a 

vector of shape design variables {V, } will not 

require the differentiation of the transformation 
matrix. For example, the following relation could 
transfer the sensitivity of the CFD load to the 
CSM grid: 

= \ T i jfV 
dv t 1 ^av; 

For the six methods described in Ref. 9, the 
transformation matrices are dependent directly 
on the shape changes. However, NURBS and 
IIM methods can be reformulated to result in 
methods with transformation matrices that are 
independent of the shape changes. Samareh 14 
proposed a reformulation for NURBS, and 
Maman and Farhat 18 proposed a general 
alternative approach for IIM. This paper provides 
reformulations for both the methods. 

The second ingredient of the proposed unified 
approach helps to reduce the number of 
transformation matrices. For a multidisciplinary 
application that involves n disciplines, the 

traditional process may require (n 2 —n) 
transformation matrices. However, some of 
these couplings are either weak or nonexistent. 
Figure 1 shows all possible interactions among 
eight disciplines; modeling all interaction 
requires 56 transformation matrices. 

The problem can be further complicated for a 
variable-fidelity multidisciplinary application. For 
example, aerodynamic loads can come from 
wind tunnel databases, or from linear 
aerodynamics, potential flow, Euler, and Navier- 
Stokes analysis codes. The following section 
provides details of an approach where the 

number of couplings is reduced from ( n 2 — n ) 
to 2 n by incorporating a CAD model. 


CAD-Based Approach 

The new element of the CAD-based approach 
was the use of a CAD geometry representation 
to reduce the number of transformations from 

( n 2 -n) to 2 n. The reduction is accomplished 
by transferring the data to a CAD geometry 
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model that serves as a common bridge or a data 
bus for sharing information among various 
disciplines, as shown in the Fig. 2. 

This approach is of obvious benefit for 
multidisciplinary applications with more than 
three disciplines. However, the method offers 
benefits even when only two disciplines are 
involved. The intermediate CAD grid isolates 
each discipline model from changes to all other 
discipline models. The overall approach has a 
strong potential for robust automation of all 
transformations, and ensures consistency 
among different renderings of the same 
configuration. 

First the data is transferred from the individual 
source discipline to the CAD model as 


Then the data is transferred from the CAD to the 
target disciplines. 


The advantage of this approach is that only 
[7 C1 ] and [r 2C ] transformation matrices have to 

be calculated: 2 n matrices instead of (n 2 — n ) . 


Transformation Matrices 


It is possible to reformulate the existing methods 
in terms of transformation matrices. This section 
presents reformulations of two existing 
alternative methods. 

NURBS-Based Interfaced 
The NURBS-based process has two steps. In 
the first step, the discipline model is mapped to 
the CAD model. Then in the second step, a 
NURBS representation is used to fit the data. 

Most CAD systems provide tools to save a CAD 
model as a NURBS representation, which then 
can be used for the mapping step. This section 
contains a brief overview of NURBS 
representation; readers should consult Refs. 14 
and 22 for more detailed discussion. A NURBS 

surface, R ( u , v) , can be represented as 


R(u,v) = -VS 

' j 

where the P tj are the locations of NURBS 
control points, the W tj coefficients are the 
weights, u and v are the parameters, 
/and /are the numbers of control points in 
u and v directions, and B ip {u ) and B jq {v) 

are the B-spline basis functions of degrees 
p and q respectively. This equation can be 

written in a compact form by combining the 
weights and basis functions into a single term as 


— ,tJ — -J— 1 

/?( M ,v)=Xc>,v)p n =[c]H 

n 

where n = i + 1 *{j - 1), and 


C„ (u,v) = C iJ (u,v) = 


W v B itP (u)B^(v) 
XlXfl,, ,(«)*, „(v) 

k 1 


Typically the data from the source discipline, 
= {/,,/ 2 ,...,/ m ,...} T . defines loads or 
deflection vectors at a discrete set of points r m . 
In the first step we project each source grid point 
r m onto the CAD NURBS surfaces and find the 
appropriate surface and the associated 
parameters, « m and v m . The projection reduces 
the number of independent variables from three, 
- *0 two, f m (u m ,v m ). The 
independent variables ( u m , v m ) are the 

parametric coordinates of the point r m on the 

NURBS surface. This information may be 
available from the grid generation process. If 
not, the grid points can be projected onto the 
original NURBS surface (see Ref. 23). 

A NURBS surface is then fitted through the data 
to form 
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/,(«,v) = [c]{f c } 

The fitted surface is based on a least-squares 
approximation (see Refs. 14-16) that minimizes 
the approximation error. The weights and the 
knot vectors of the CAD NURBS surface can be 
used, or the user can specify weights and knot 
vectors. The degrees of NURBS approximation, 
p and q , and the knot vector distribution affect 

the smoothness of the least-squares 
representation. The minimization error E can 
be written as 

}- [cite }} r (te }- [cite }} 

The least-squares form of the above equation 
can be expressed in a matrix form as 

te}=te,ite} 

where 

tei= [carter 

A generic high-speed civil transport (HSCT) 
geometry was used to demonstrate the 
algorithm 14 . This geometry was made of three 
surfaces: fuselage, inboard wing, and outboard 
wing. Figure 3 shows the original (undeflected) 
NURBS surfaces, the deflected CSM grid, and 
the deflected NURBS surfaces. To test the 
limits of the deflection-transfer algorithm, the 
CSM grid had a large and unrealistic deflection. 

Because the CSM grid is generally coarser than 
the CFD surface grid, the interpolation approach 
for deflection transfer may produce a 
discontinuous CFD surface grid. However, the 
NURBS-based approach maintained the 
smoothness of the geometry for deflection 
transfer. This is a major advantage for the 
NURBS based approach. On the other hand, for 
transferring the integrated quantities, the 
process in its present form did not guarantee 
conservation of forces and moments. Further 
studies need to be performed to determine 
whether the advantages from the smoothness of 
the NURBS approach outweigh the lack of 
conservation. 


[ nterpolation-Based Algorithm 
We followed the algorithms proposed in Refs. 18 
and 21. The first step is to map each source and 
target grid point to the CAD model as described 
earlier, and map the CAD model to the source 
and target disciplines. We used a discrete 
representation of the CAD model, which is 
readily available from most CAD systems. Each 
source and target point was mapped to the CAD 
model, and the parametric coordinates (u, v) of 
the source and target grid points were 
determined from the mapping process. The 
second step is to transfer the data from the 
source discipline to the target discipline. In this 
step, the parametric coordinates were used to 
transfer the data to the CAD model: 

i 

where the terrniV, was the FE shape function. 

This equation was applied for each source grid 
point. The resulting set of equations was 
assembled into 

te)=te,ite} 

The elements of the transformation matrix were 
made of the FE shape functions. A similar 
process was used to transfer the data from the 
CAD model to the target discipline. As 
demonstrated by other researchers, 7, 8 the IBA 
was very effective in transferring the scalar and 
vector quantities. 

The original 18 and the modified IBA 21 were used 
to transfer the aerodynamics data for a blended- 
wing body and an HSCT model. Figure 4 shows 
the result of transferring sensitivity derivative 
data from a CFD grid to a CSM grid for an HSCT 
model using the original algorithm. The data is 
the sensitivity derivative of pressure with respect 
to the leading-edge sweep angle. 

Figure 5 shows the result of transferring the 
pressure distribution for a blended-wing body 
using the original algorithm. The figure shows 
the pressure contours on the CFD grid and the 
pressure contours transferred from the CFD grid 
to the CAD model. It also shows the pressure 
contours transferred from the CAD model to the 
CSM grid. As expected, the original IBA 18 did not 
maintained the conservation of forces and 
moments. The error for the integrated forces 
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was approximately 5%. This error would be 
higher when the models have dissimilar levels of 
geometry details. 

For transfer of the integrated quantities, the 
modified IBA had a slightly different 

implementation (see Ref. 21), which conserves 
forces and moments. Figure 6 shows the load 
vectors on the CFD grid and the load vectors 
transferred from the CFD grid to the CAD model. 
It also shows the load vectors transferred from 
the CAD model to the CSM grid. As expected 21 , 
the transfer process maintained the 

conservation of forces and moments within the 
machine’s accuracy. This process can also be 
used to transfer deflections from a CSM grid to a 
CFD grid. 

Figure 7 shows the result of transferring the 
aeroelastic deflections from a CSM grid to a 
CAD model and then to the CFD grid. The top- 
left figure shows the deflected CSM grid. The 
deflection is transferred from the CSM grid to the 
CAD model, as shown in Fig. 7 (middle-left 
figure). Then the deflection is transferred from 
the CAD model to the CFD grid (bottom-left 
figure). The right portion of the figure shows the 
deflected CSM and CFD grids, and they are 
right on top of each other. A closer look at the 
Fig. 7 demonstrates the benefit of the three- 



g eundpsia 


Fig. 1 Multidisciplinary interactions. 


dimensional data transfer. The winglet is not 
merely shifted up — it has followed the CSM 
deflection without neglecting the x and y 
components of the deflection. The traditional 
data transfer will only include the vertical 
displacement ( z coordinates). 

Summary 

We have presented a unified approach for 
transferring information for a multidisciplinary 
application with n disciplines. This method has 
two essential ingredients. First, the data transfer 
process between two disciplines was modeled 
by a transformation matrix. Second, the CAD 
model with consistent geometry was used to 
reduce the number of interactions from 

(n 2 -n) to 2 n. This unified approach was 
specially tailored for application environments 
where the geometry is created and managed 
through a CAD system. Results were presented 
for a high-speed civil transport and a blended- 
wing body. The approach provides a framework 
where data transfer among disciplines can be 
accomplished consistently, and potentially with a 
high degree of automation. Comparisons 
between the NURBS and IBA using CAD and 
transformation matrices are continuing. 



Fig. 2 Unified multidisciplinary interactions. 
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Fig. 3 Aeroelastic deflection transferred from CSM grid to a CAD model. 
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Fig. 7 Transfer of aeroelastic deflection from CSM grid CFD grid. 
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